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Abstract. 

We revise the Lewis-Riesenfeld invariant method for solving the quantum time- 
dependent harmonic oscillator in light of the Quantum Arnold Transformation 
previously introduced and its recent generalization to the Quantum Arnold-Ermakov- 
Pinney Transformation. We prove that both methods are equivalent and show the 
advantages of the Quantum Arnold-Ermakov-Pinney transformation over the Lewis- 
Riesenfeld invariant method. We show that, in the quantum time-dependent and 
damped harmonic oscillator, the invariant proposed by Dodonov & Man’ko is more 
suitable and provide some examples to illustrate it, focusing on the damped case. 


1. Introduction 

The Lewis-Riesenfeld invariant method [1, 2] is a technique that allows to obtain a 
complete set of solutions of the Schrodinger equation for a time-dependent harmonic 
oscillator in terms of the eigenstates of a quadratic invariant. This quadratic invariant, 
the Lewis invariant, is built using an auxiliary variable that satisfies the Ermakov 
equation [3, 4, 5, 6]. 

The Quantum Arnold-Ermakov-Pinney Transformation (QAEPT) [7] is a unitary 
transformation that maps solutions of a Generalized Caldirola-Kanai [8, 9, 10] 
Schrodinger equation into solutions of another Generalized Caldirola-Kanai Schrodinger 
equation. In particular, one of the systems can be the standard harmonic oscillator and 
the other a time-dependent harmonic oscillator, and in this case we shall show that the 
Lewis-Riesenfeld invariant method is recovered. 

The idea of using invariants to solve equations is rather old, going back to S. Lie 
(1883) [11] who showed that a second order differential equation has the maximal group 
of symmetries if the differential equation is up to third order in the derivative, and the 
coefficients satisfy certain relations [12, 13]. 

V.P. Ermakov (1880) [3, 4] showed that the general solution of the non-linear 
equation 

b + c 


( 1 ) 
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where ui a is an arbitrary constant, can be obtained from two independent solutions y \, y 2 
of the corresponding linear equation: 

ij + uj 2 (t)y = 0 (2) 


by: 

b 2 = ci y\ + c 2 yl + 2c 3 y 1 y 2 , (3) 

with C 1 C 2 — c 3 = c <Jq. Similar results were derived independently by W.E. Milne (1930) 
[5] and E. Pinney (1950) [6]. 

H.R. Lewis (1967) [1] obtained a classical and quantum quadratic invariant for a 
time-dependent harmonic oscillator, of the form: 

1 • 1 

x L = —(bp-m bx ) 2 + -rmj 2 —, (4) 

where, again, ojq is an arbitrary constant with dimension of frequency, in terms of an 
auxiliary dimensionless function b(t) satisfying the Ermakov equation (1). Note that 
classically x satisfies the equation of motion 

x + uj 2 {t)x = 0 . (5) 

The pair of equations (1) and (5) is denoted an Ermakov system. They are uncoupled 
(given l a(t), they can be solved independently for x and b ), although it has been 
generalized to coupled equations and to higher dimensions [14, 15]. 

The reader might wonder about the comparison of b (dimensionless), satisfying (1) 
and providing the invariant (4), and a function p satisfying the usual Ermakov-Pinney 
equation p + ou(t) 2 p = which has the dimensions of the square root of time (see e.g. 
[1]). The relation between b and p is simply b = y/uJop. Also, the relation between Z L 
and the invariant / in [1] is simply Zl = (jJqI, i.e. Zl has dimension of energy meanwhile 
/ has dimension of action. The arbitrariness in the choice of b and Zl was already noted 
in [1, 2] and we make it explicit for later convenience (see also [16], where the authors 
use the same convention). 

Lewis & Riesenfeld (1969) [2] used the eigenvectors of the quantum version of this 
quadratic invariant Zl written in terms of the auxiliary function b{t) satisfying the 
Ermakov equation to obtain solutions of the Schrodinger equation for a time-dependent 
harmonic oscillator. For this purpose, an extra time-dependent phase e* ^ had to be 
added to the eigenfunctions in order to satisfy the Schrodinger equation. They did not 
considered damping (or time dependent mass) and they supposed that the quadratic 
invariant has discrete spectrum. 

V.I. Arnold (1978) [17], in the context of symmetries of second order ordinary 
differential equations, introduced the term straightening for the linearization studied by 
S. Lie and considered the case of Linear Second Order Differential Equations (LSODE): 

x + fx + u 2 x = A , (6) 
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where /, uj and A are time-dependent functions, giving explicitly the transformation for 
this case: 

A : RxT -x IxT 
(x,t) (k,t) 


(V 


with 


r = 


Ul(t) 

u 2 (t ) ’ 


K = 


X — U p (t ) 

u 2 (t) 


( 8 ) 


where T and T are, in general, open intervals; U\ and u 2 are independent solutions of 
the homogeneous LSODE, u p is a particular solution of the inhomogeneous LSODE, 
and W(t) = u\u 2 — u\u 2 = is the Wronskian of the two solutions]:. 

Under this transformation, the classical equation of motion (6) transforms as: 

A W 


x + fx + to x = A 


K = 0 . 


U n 


( 9 ) 


Thus, the Arnold transformation maps patches of solutions of the LSODE system into 
patches of free particle trajectories. 

For convenience, we shall impose the canonicity conditions (see [10]): 

wi(0) = u 2 (0) = w p (0) = Mp(0) = 0 , fti(0) = w 2 (0) = 1 • (10) 


These conditions play a crucial role in the physical interpretation of quantities mapped 
from one system into the other through the Arnold transformation. More precisely, if 
k(t) and n (r) are the conserved position and momentum for the free particle (verifying 
that k(0) = k and 7r(0) = n = mk), then the transformed quantities through the Arnold 
transformation are the conserved position x(t) and momentum p(t) in the LSODE 
system (verifying x(0) = x and p(0) = p = mx). 

Dodonov & Man’ko (1979) [18, 19] (and Malkin-Man’ko-Trifonov [20] (1969) 
without considering damping) computed the coherent states for the Generalized 
Caldirola-Kanai model (the quantum version of a general LSODE system), whose 
Hamiltonian is 

V 2 1 

Hgck = + ( -mu 2 x 2 - mAx)e f , (11) 

2m 2 

using first-order invariants as annihilation and creation operators. The number operator 
associated with these annihilation and creation operators is a quadratic invariant that 
will be denoted the Dodonov-Man’ko invariant Xdm- Later, other authors have used 
first-order invariants to solve time-dependent problems [21]. 

Hartley & Ray (1981) [22] and Lewis & Leach (1982) [23] generalized the 
construction of the Lewis invariant to some non-linear systems. 

Pedrosa (1987) [24] constructed the Lewis invariant for the Ermakov equation with 
a damping term using canonical transformations. 


f If x represents position and t represents time, then u\ has dimension of time, u 2 is dimensionless and 
Up has dimension of length. 
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V. Aldaya et al. (2011) [10] extended to the quantum case the Arnold 
transformation and denoted it the Quantum Arnold Transformation (QAT): 

A : l-l t —» U° 

(fi(x,t) 1 —> <p(k,t) = A (4>(x, t)) 

“ 2(0 -2 


( ,- i m 1 “2 0) -2 \ 

sjU2(t) e 2 h WO) “ 2(0 (j)(x,t)j 


( 12 ) 


ffere A* is defined as A*(f(x,t )) = /(A _1 (fi;, r)), TR is the Hilbert space of solutions of 
the Generalized Caldirola-Kanai Schrodinger equation at time t, and is the Hilbert 
space of solutions of the Schrodinger equation for the free Galilean particle at time 
r, where t and r are related by the Arnold Transformation. Note that the QAT 
transforms solutions of the time-dependent Schrodinger equation of the Generalized 
Caldirola-Kanai system into free-particle wave functions, and that this is achieved by 
applying the Arnold transformation together with multiplying the wave function by a 
suitable phase and rescaling factor. These factors also renders the QAT unitary [10, 7]. 

Some applications of the QAT were given in [25], where states from the harmonic 
oscillator were mapped into the free particle giving rise to Hermite-Gauss and Laguerre- 
Gauss wave packets; in [26, 27], where processes of Release and Recapture of a particle by 
a harmonic trap were studied using the QAT; and in [28], where the QAT, which is a local 
diffeomorphism in time, is extended beyond the “focal” points, correctly reproducing 
the change in phase of the wave function (Maslov correction, see for instance [29]). 

Castanos, Schuch & Rosas-Ortiz (2013) [30] constructed coherent states for 
different models (time-dependent and non-linear Hamiltonians) through complex Riccati 
equations and found the corresponding Lewis invariants. 

Since its introduction, the Lewis invariant and its associated Ermakov equation 
entered an inflationary scenario with applications in many areas. One of the most 
remarkable ones are the applications in Bose-Einstein Condensates (BEC) [31, 32], 
where a transformation similar to that of Arnold (and known as scaling transformation in 
this context) taking the time-dependent harmonic trap in the Gross-Pitaevskii equation 
into a stationary one is applied. Although they did not use the Lewis invariant, the 
scaling parameter satisfies the Ermakov equation (1). In this context, eq. (1) also 
appears in [33]. 

Recently, the Lewis-Riesenfeld invariant method has been used to inverse engineer 
short-cuts to adiabaticity [16], to speed up cooling processes and transport in 
electromagnetic traps and BECs, and to manipulate states in wave-guides [34], where 
the relation with Generalized Caldirola-Kanai systems has been established [35]. The 
main idea here is to design a Lewis invariant satisfying the property of commuting with 
the Hamiltonian at initial and final times, and this can be achieved by building up a 
function b satisfying certain boundary conditions and then determining, through the 
Ermakov equation, the time-dependent frequency that should be applied in order to 
take the system from the initial state to the desired final state without affecting to the 
population of the levels. 
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Another recent application of the Lewis-Riesenfeld invariant method is in 
mesoscopic RLC electric circuits [36], where the quantum evolution (even in the case of 
time-dependent R(t), L(t ) and C(t) and source term) is described. 

The content of the paper is as follows. In Sec. 2 we revise the Lewis-Riesenfeld 
method and explain it in terms of the QAT and the QAEPT, showing that the use of 
the Dodonov-Man’ko invariant is more appropriate for damped systems. In Sec. 3 the 
examples of the Caldirola-Kanai and the Hermite oscillators are studied in detail. 


2. The Lewis-Riesenfeld method in light of the Quantum Arnold 
Transformation 


In their original paper Lewis & Riesenfeld [2] provided a method to obtain a family 
of exact wave functions for the time-dependent harmonic oscillator spanning the whole 
Hilbert space. In a first step the method looks for an invariant, Hermitian operator Zl, 
a task which can follow the lines of [1]. Imposing the invariance condition§ 


dZL 

di 


f AP( ( )A = o. 


(13) 


where H(t) is the Hamiltonian for the time-dependent harmonic oscillator and assuming 
the most general quadratic invariant, they arrived at the quantum version of (4), where 
the auxiliary function b{t) satisfies the Ermakov equation (1). In this equation oo 0 is just 
an arbitrary constant. The possibility exists of giving a generic form for the invariant 
(i.e. quadratic or linear, Hermitian or complex combinations of basic operators p and 
x, etc.) and then solving for the coefficients to fulfill (13). 

The second step in the method is realizing that finding eigenfunctions (j) s (x, t) of Zl, 
Zl 4> s (x,t) = X s (j) s (x,t), amounts to finding solutions of the Schrodinger equation except 
for a time-dependent phase, which must be computed. That is, solutions il) s {x,t) of the 
Schrodinger equation 



h 2 d 2 i/j s 
2 m dx 2 


+ -muj 2 (t)x 2 'ip s 




(14) 


may be of the form 


t) = e iQs(i V s (:E,l), 


(15) 


where a s (t) satisfies 

- H(t)j </> a (x, t ). (16) 

§ Although we are working in the Schodinger picture of Quantum Mechanics, i.e. wavefunctions depend 
explicitly on time and common operators like position x and momentum p do not depend on time, other 
quantum operators may depend explicitly on time. This is precisely the case of the Hamiltonian for 
non-conservative systems and in general for invariant operators. 
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That is just a nice consequence of the fact that Z L applied on a solution of (14) is 
again a solution: 



(17) 


The phase can be solved in terms of h to give: 



(18) 


Two observations can be made. First, the way in which the eigenfunctions (j> s (x,t) 
of the invariant are found are left to the ability of the user of the method. In this 
respect, some authors have developed a unitary transformation from the Hamiltonian of 
the simple harmonic oscillator into the invariant [37] (resembling very much the QAEPT, 
see below). For this purpose, transforming the time-independent Schrodinger equation of 
the simple harmonic oscillator, including the wave functions, would do the trick. Second, 
although the Lewis-Riesenfeld method can provide all (quadratic) invariants for the time 
dependent harmonic oscillator (by taking all possible solutions of the Ermakov equation) 
it does not provide insight on their physical interpretation (like their spectra). In our 
case, we provide a method that allows a neat physical interpretation of the invariants 
since each one preserves the same character as in the harmonic oscillator. 

In that sense the QAT and its generalization the QAEPT turn out to be very useful. 

2.1. The Quantum Arnold-Ermakov-Pinney transformation 

The QAEPT is obtained when two different LSODE-systems are related by QATs with 
the free-particle system as an intermediary, that is, when a QAT and an inverse QAT 
are composed. That was shown in [7]. As the QAT, the QAEPT relies on the symmetry 
structure of the systems of the Generalized Caldirola-Kanai type (see [10]). 

In the QAEPT, it is the full set of invariant operators and the corresponding 
eigenstates including the time dependence (with no need of searching for a phase as 
it is given by the transformation) which is mapped from the simple harmonic oscillator 
system into a Generalized Caldirola-Kanai system. The interpretation of the eigenstates 
may be the same on both sides of the mapping. 

Let Ai and A 2 denote the Arnold transformations relating the LSODE-system 1 
and LSODE-system 2 to the free particle, respectively, then E = Af 1 A 2 relates LSODE- 
system 2 to LSODE-system 1. E can be written as: 


E : M x T 2 —y M x Tj 

(x 2t t 2 ) (xi, ti) = E(x 2> t 2 ). 


(19) 


The explicit form of the transformation can be easily computed by composing the 
two Arnold transformations, resulting in: 
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where b(t 2 ) 


( 2 ) 

u f satisfies the non-linear SODE: 
“2 hi) 


b + f‘2 b + u 2 b 



( 21 ) 


and where refers to the 7 -th particular solution for system 7 : Wj, fj and c Oj stand for 
the Wronskian and the LSODE coefficients for system 7 ; and the dot means derivation 
with respect to the corresponding time variable. If all uf i satisfy the canonicity 
conditions ( 10 ) then 6 (i 2 ) satisfies the corresponding canonicity conditions 


6 ( 0 ) = 1 , 6 ( 0 ) = 0 . 


( 22 ) 


Equation (21) constitutes a generalization of the Ermakov equation. That equation, 
together with the LSODE of system 2, is a generalized Ermakov pair [14, 15]. Also, any 
(quadratic) conserved quantity, which is shared by the two LSODE-systems, constitutes 
a generalized Lewis invariant. Equation (21) actually defines a generalized Arnold 
transformation, to be named (classical) Arnold-Ermakov-Pinney transformation , which 
transforms solutions of the LSODE 1 into solutions of the LSODE 2. 

The quantum version of the Arnold-Ermakov-Pinney transformation, E, can be 
obtained by computing the composition of a QAT and an inverse QAT to give: 


E : 




fi(x 2 ,t 2 ) 


<p(xi,ti) = E (0(x 2 ,t 2 )) 


= E* 



1 rn 1 b{t 2 ) 9 

2 h w 2 (t 2 ) b(t 2 ) 2 0(x 2 ,£ 2 ) 


(23) 


The Quantum Arnold-Ermakov-Pinney transformation (QAEPT) maps solutions 
of a Generalized Caldirola-Kanai Schrodinger equation into solutions of a different, 
auxiliary Generalized Caldirola-Kanai Schrodinger equation, and by construction it is 
also a unitary transformation. The auxiliary system might be, in particular, the one 
corresponding to a harmonic oscillator with frequency cui(t) = ojq and f± = 0. In this 
case eqs. ( 20 ) and ( 21 ) reduce to: 


and 


X2 

Xl ~ b(t 2 ) 



b(ty- ’ 


6 + f 2 b + 6^2 b 


W 2 2 ^ 

b 3 


(24) 

(25) 


2.2. Ermakov System and interpretation of the Lewis Invariant 


Consider the particular case where LSODE-system 1 is a harmonic oscillator (cui(ti) = 
oj 0 and f\ = 0), which can be described by the Hamiltonian 


Hho — 


Pj 

2m 


2 2 


+ -mujnX 




(26) 


and LSODE-system 2 is a time-dependent harmonic oscillator with frequency o; 2 (f 2 ) = 
oj{t ) and / 2 = 0, with Hamiltonian H(t ) given by (14). Then, expressions (21) and 
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(25) simplify to (1). Obviously, for cuo = 0 the Arnold-Ermakov-Pinney transformation 
reduces to the ordinary Arnold transformation, i.e. E = A. 

Now, note that LSODE 1 Hamiltonian, Hho , is conserved, and that it is so on both 
sides of the transformation E , given by (see (20)): 

X f 1 1 

x\ — — , t\ — I —dt = — arctancu 0 T , (27) 

b J b 2 u 0 

where r denotes the (common) time in the free particle given by the Arnold 

transformations A\ and A 2 . Also, b(t) = 1 + uJqT 2 = \J ( u^) 2 + satishes 

the Ermakov equation (1) together with the canonicity conditions (22). 

It should be stressed that b{t) never vanishes, otherwise the Wronskian W 2 {t) of 
the two independent solutions would also vanish. And since in the quantum case the 
time t\ appears in the form e~ lUotl = e “ !arctan (‘ doT ) ) this expression is well-defined for 
all times (even in the case where r has singularities). This means that the QAEPT 
transformation is well defined for all times (i.e. £ : 1 x K ->1x1). This is an 
important advantage with respect to the QAT, that was defined only locally in time. 

Computing the momentum pi = mx\ = = m(xb — bx ), we can 

write Hho hi variables corresponding to system 2: 

H H o = ^{pb-mbx)‘ 2 + ^mw 0 2 (^) 2 = X L . (28) 

That is easily recognized as the usual Lewis invariant Zl. Thus, we have found a way to 
characterize it through the Arnold-Ermakov-Pinney Transformation: Zl corresponds 
to the conserved quantity Hho imported from the simple harmonic oscillator, which is 
used as an auxiliary system. Because the auxiliary system is arbitrary, Zl is conserved 
for any Cti 0 , provided (1) is satished. Note that, in order to establish the identification 
Hho = 27 it is essential to impose the canonicity conditions (22). 

Using the explicit form of the inverse E~ l of (23) in this case, it is straightforward to 
arrive at solutions </>(x, t) of the Schrodinger equation of the time-dependent harmonic 
oscillator in terms of solutions of the Schrodinger equation for the simple harmonic 
oscillator (p(xi,ti): 

= ’ ( 29 ) 

where b is any solution of (1) satisfying the canonicity conditions (22). Note 
that, if ip(x i,U) is chosen to be, for instance, an eigenfunction of the quantum 
operator corresponding to (26), HhOi then the transformed wave function (p(x,t) is 
an eigenfunction of the quantum operator Zl corresponding to the invariant (28) (the 
explicit form of such operators is easily obtained from their classical counterpart by the 
canonical quantization prescription). That shows that Zl has discrete spectrum. 

In the wave functions (29), two phases can be distinguished: the one corresponding 
to the transformation itself, explicit in (29), and the phase mapped from 
(which is the only time dependence for stationary states in the harmonic oscillator) into 
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e *( n+ 2 * , 0)2 _ The latter accounts for the phase of the Lewis-Riesenfeld method. The 
former accounts for the phase (and the factor) which appears, for instance, in [37]. 

Regarding the canonicity conditions (22), they play an important role in short-cuts 
to adiabaticity processes for time-dependent harmonic oscillators, see [16], since they 
imply that Xl commute with the Hamiltonian at the initial time t = 0. If we further 
impose b = 0, then oj 0 = 0 ^( 0 ) holds, and the Hamiltonian at t — 0 will coincide with 
the invariant Xl at t = 0. In the following, we shall assume that the Lewis invariant JR 
verifies these conditions. 

The same process can be repeated for any other operator representing an invariant 
in the simple harmonic oscillator (LSODE 1), showing the usefulness of the QAEPT to 
perform quick computations. 

In conclusion, it is the full set of invariant operators and the corresponding 
eigenstates (with no need of searching a phase) what is mapped from the simple system 
into the Generalized Caldirola-Kanai system through the QAEPT. Also, a word of 
caution is in order: the Hamiltonian operator of one system is not mapped into the 
Hamiltonian operator of the other system, which may not be invariant itself. 


2.3. The Lewis-Riesenfeld (Dodonov-Man’ko) invariant method for the Generalized 
Caldirola-Kanai oscillator through the QAEPT 


In a Generalized Caldirola-Kanai system, the easiest way to find eigenstates of an 
invariant operator and its eigenfunctions as solutions of the Generalized Caldirola-Kanai 
Schrodinger equation, is to focus on an auxiliary system (the harmonic oscillator in the 
previous subsection) with its Hamiltonian being the invariant operator and perform 
the QATs or QAEPT necessary to map the Schrodinger equation of such auxiliary 
system into the Generalized Caldirola-Kanai Schrodinger equation. In this process, 
Hamiltonian is not mapped into Hamiltonian, but conserved operators into conserved 
operators are. That procedure takes advantage of the fact that the eigenstates of the 
harmonic oscillator Hamiltonian have a very simple time dependence, as has just been 
noticed in the previous Subsection. 

Let us now describe a different way of constructing an invariant. The idea is 
to consider any linear combination of quadratic invariants in such a way that its 
eigenfunctions solve the Generalized Caldirola-Kanai Schrodinger equation. The most 
general invariant can be written in the form [ 10 ]: 


X 


1 1 9 a „ 7 XP + PX 

— P 2 + - mij 2 X 2 + ---- 

2m 2 2 2 


(30) 


where u) and 7 are arbitrary real numbers and X, P are conserved position and 
momentum operators satisfying that at t — 0 coincide with the usual x,p , namely 



On the Lewis-Riesenfeld (Dodonov-Man’ko) invariant method 


10 


The new invariant X plays now the role of Hho i n this more general setting. 
The eigenfunctions of this operator, solutions of the Generalized Caldirola-Kanai 
Schrodinger equation, are [10]: 





U2 — iLlu\ 


2 Tyrmx 

e 


Cl 


^2 ) 
u 2 W J 


C\D U 



+ C2D-1- 



(32) 


where C\ and C2 are arbitrary constants, Cl = y / Cj 2 — ^ and v is in general a complex 
number, D v are the parabolic cylinder functions and T(z) is the Gamma function [39], 
Mi and w 2 are solutions of the LSODE equation corresponding to the given Generalized 
Caldirola-Kanai oscillator, w 2 = w 2 — 714/2 and the function 

6 ( 2 ) = \J u\ + Cl 2 u\ = \J (w 2 — / fUi/2) 2 + Cl 2 u\ (33) 

plays the role of 6 ( 2 ), i.e. satisfies the Generalized Ermakov equation (25) with jjq = LI, 
but with different initial conditions, namely 5(0) = 1 and 6 ( 0 ) = — 

As in the case without damping, neither h{t) nor b{t) vanish, otherwise the Wroskian 
W ( t ) of the two solutions u\ and m 2 would also vanish. 

Thus, the invariant X can be seen as Hho mapped from the harmonic oscillator with 
frequency 12 through a QAEPT (24) characterized by b[t) satisfying (25) with uq = 12. 

The associated spectrum of X is 

\ v = hCl (y + —). (34) 

The integer, real or complex character of v depends on the value of 12 (and this in 
turns depends on the particular values of Cj and 7 ). See [18] for a discussion in the case 
of the damped harmonic oscillator. 

The choice Co = a; 2 (0) and 7 = 0 leads to a generalized Lewis invariant X L where 
the arbitrary frequency has been chosen as cuo = 1^2 ( 0 ), providing an invariant that 
commutes with the Generalized Caldirola-Kanai Hamiltonian (11), and in fact coincides 
with it, at t = 0. This invariant could be useful in short-cuts to adiabaticity processes 
for damped systems or with time-dependent mass (like in waveguides [35]), where an 
invariant commuting with the Hamiltonian is needed. 

For the damped harmonic oscillator with constant co and 7 (see Sec. 3.1), there 
exists a different choice Cj = iv and 7 = 7 leading to the only quadratic invariant 
X = Ldm, the Dodonov-Man’ko Invariant, whose unique, explicit time dependence is 
through the Wronskian W(t) (like the Caldirola-Kanai Hamiltonian itself). Even more, 
the coherent states associated with this invariant (through a factorization of the form 
Zdm — + AA^), A being conserved creation-annihilation operators) are the 

only ones with minimal, time-independent uncertainty relations (see Dodonov & Man’ko 
[18]). 
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In conclusion, when studying Generalized Caldirola-Kani systems with the Lewis- 
Riesenfeld invariant method or with the more general method given by the QAEPT, a 
different choice for the invariant operator than Zl should be made, in particular Z (with 
suitable coefficients Co and 7 ) may be more appropriate (like the case of Z DM for the 
damped harmonic oscillator). We shall denote this invariant a Generalized Dodonov- 
Man’ko invariant Zgdm- 


2-4- Engineering a suitable QAEPT to build a Generalized Dodonov-Man’ko invariant 

Once the general setting has been established, let us apply the method in a suitable 
way to obtain a proper invariant Zgdm- The previous discussion on the choice of 
an appropriate invariant, together with the corresponding analysis of the Caldirola- 
Kanai oscillator (with constant damping and frequency, see below), suggest that it 
would be helpful to construct a QAEPT from a Generalized Caldirola-Kanai system 
2 to a yet undetermined Generalized Caldirola-Kanai system 1 , but satisfying certain 
requirements, implemented in the choice of b(t 2 ). In other words: we look for an auxiliary 
system to help solving a Generalized Caldirola-Kanai oscillator in such a way that the 
QAEPT is as simple as possible. 

In particular, it is easy to check that choosing b = W 2 ' , the Generalized Caldirola- 
Kanai system 1 is a time-dependent harmonic oscillator, i.e. it has no damping term. 

In fact, with b = W 2 the corresponding QAEPT is given by: 


1 1 = 


W2 

b 2 


dt 2 = t 2 = t , 


x 2 

Xi = — = 


X 2 


w. 


1/2 ’ 


and 


( t>(x 2 ,t ) = W 2 ' e 


-2 , x 2 

<^(777173.0 - 


w. 


(35) 


(36) 


The wave function ip satisfies the Schrodinger equation for a time-dependent harmonic 
oscillator with frequency 

2W 2 W 2 - 3W 2 = ^ (t)2 l,o 1 , 

4" 


uq(f) 2 = u 2 {tf + = u 2 {tf - t / 2 2 - y 2 , 


4W 2 


(37) 


and the Wronskian for system 1 turns to be W\ = 1, i.e. the auxiliary system 1 
is not damped. However, 07 (t) is not arbitrary but specifically designed to simplify 
the mapping. Therefore, it is straightforward to map results (such as invariants) and 
computations from the known, auxiliary system 1 to the Generalized Caldirola-Kanai 
system 2 . 

In this case the canonicity conditions (22) are not satisfied in general, since b( 0) = 1 
but 6(0) = \W{ 0) (see the examples on Sec. 3). 

Note that transformation (35), together with its corresponding extension to velocity, 
is nothing other than a generalized version [24] of the time-dependent canonical 
transformation that removes the damping in the damped harmonic oscillator [46] (see 
also [47, 48]). 





On the Lewis-Riesenfeld (Dodonov-Man’ko) invariant method 


12 


Constructing the Lewis invariant Xl for this time-dependent harmonic oscillator, 
and mapping it back to our Generalized Caldirola-Kanai system 2 through the previous 
QAEPT, leads to a Generalized Dodonov-Man’ko invariant Xgdm for system 2 . More 
precisely, the invariant can be written as: 


Xqdm — 


P 2 
2m 


-m 



5/2(0)) X 2 + 


M 0 ) XP + PX 
2 2 


(38) 


We shall provide examples of this construction in the next section. 


3. Examples 

Let us discuss some examples of damped systems, where the previous ideas can be 
applied to construct invariants and simplify analytical computations. 

3.1. Caldirola-Kanai oscillator 

The simplest example that can be studied is the quantum damped harmonic oscillator, 
also known as Caldirola-Kanai oscillator [ 8 , 9]. From the point of view of invariant 
operators, it was first studied by Dodonov & Man’ko [18], who constructed first-order 

invariants in the form of conserved creation and annihilation operators and derived 

a basis of number states and a family of coherent states satisfying minimal, time- 
independent uncertainty relations. The classical equation of motion is given by (we 
shall not consider the external force term): 

x + qx + u 2 x = 0 , (39) 

where oj and 7 are constants. Even though this is a linear equation with constant 
coefficients, the system is not conservative since the Hamiltonian describing this equation 
is time-dependent: 

Hck = e' 7 * W . (40) 

There is an old controversy with the quantum version of the Caldirola-Kanai 

oscillator concerning the dissipative character of this system [40]. The main drawback 
is that the evolution is unitary for all times and there is no loss of coherence, something 
that it is considered inherent to a quantum dissipative system. Some proposals have 
been made to address that paradoxical situation [41] (see also [42]). 

An alternative physical interpretation to the damping term in (39) is that of a 
time-dependent mass, that is, the mass is actually of the form: 

1 

m(t) = me 1 * =► H C k = + ^rn(t)uj 2 x 2 . (41) 

Thus the Caldirola-Kanai oscillator describes an oscillator whose mass is growing 
exponentially. 
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The solutions of the classical equations (39), satisfying the canonical conditions 
Wi(0) = 7 / 2 ( 0 ) = 0, 7 / 2 ( 0 ) = tti(0) = 1 (see [10]) are: 


x , sin Lit 


ui(t) = e 2 * 


n 


U 2 (t) — e 2 1 (cos Lit + ^ sin f Itj , 


(42) 


where 'uj 2 — ^ and the Wronskian is given by W(t) = t/i(1)t/2(1) — U\{t)u 2 {t) = 

e _7f = m/m(t). The solutions are oscillatory if u > ^ (underdamping) and have a 
single zero for u < 7 (overdamping) and cu = 7 ) (critical damping). In this last case the 
solutions are u\{t) = le~ and 7 / 2 ( 1 ) = e'a* (l + t[1). 

Let us restrict to the underdamping case. From the previous solutions the Arnold 
transformation (7) is given by: 

sin Qt 


r = 


7/1(1) 

7/2(1) 

X 


COS Lit + si 11 ^1 


X 


(43) 


hi = 


7 / 2 ( 1 ) e (cos 111 + 2 ^ sin 111) 

The Arnold-Ermakov-Pinney transformation mapping the Caldirola-Kanai 
oscillator into the standard harmonic oscillator (and its corresponding quantum version) 
is also easily derived, resulting in: 


1 ' = 


x = 


Wit ) , 1 

dt = — arctan ujqt 
b(ty uj 0 


x 


(44) 


Kt)’ 


where ojq is the (arbitrary) frequency of the auxiliary harmonic oscillator, and 


6(1) = U 2 {t)\Jl + C^r 2 = \Ju 2 {t) 2 + LU$Ui(ty 


-Z+ 
= e 2 


cos Lit + ^ sin Lit 


, sin" Lit 


(Wn 


LI 2 


(45) 


Note that 6(1) never vanishes since this would imply u\ (1) = 7 / 2 ( 1 ) = 0 at some time 
instant, and thus the Wronskian would also be zero at that time instant, contradicting 
the fact that u± and 7/2 are independent solutions of (39). Note also that 6(1) satisfies 
the canonicity conditions (22). 

Let us construct an invariant for the Caldirola-Kanai oscillator. The first possibility 
is to construct a generalized Lewis invariant II with u = 10 and 7 = 0: 


X L = P 2 + ]-muj 2 X 2 , 

2m 2 


(46) 


leading to the invariant constructed by Pedrosa [24]. 

According to the discussion in Sec. 2 , for a damped system it is more appropriate 
to use the Dodonov-Man’ko invariant, which is built using the solution 7/1 (l) and the 
(noncanonical) solution 7/2(1), with Cj = u and 7 = 7 : 


1 . sin Lit 


7 / 1 ( 1 ) — e 2 * 


LI 


W 7 , 

7/2(1) = 7/2(1) — — 7 /i(l) = e~2* cosDl. 


(47) 
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Therefore h{t) = ^/ft 2 ( t ) 2 + ^ 2 Mi(t ) 2 = e a 4 , which also never vanishes. The 
Dodonov-Man’ko invariant (30) for these solutions is given by: 


X DM — x —^ T + ~ 


7 XP + PX 


_ _ _ _ , — H(jx -\— 

2m 2 2 2 2 2 

It can be shown that, for the specific choice u; = cu and 7 = 7 made above, the 
invariant Tdm coincides with that provided by Nassar [44, 45]. It also coincides with 
the invariant discussed by Cervero & Villarroel [38]. 

The eigenstates of this invariant, solutions of the Caldirola-Kanai Schrodinger 
equation (note that the time-dependent phase appearing in the Lewis-Riesenfeld 
invariant method is included in the QAEPT), are derived from (32): 


7 xp + px 


(48) 


< K{x,t ) = 


a/ \Z2tt2 7 


D —i (n!l+ i (n+i 2)) t — H (0+ i % )e yt x 2 


Hr, 


n\ 


mOt 

~h 


e'Xx 


,(49) 


with eigenvalues Kl{n+ |), with n — 0,1,... (assuming D > 0). Thus, the result of [18] 
is recovered. 

Note that in deriving (49) it has been crucial to choose the solutions Ui(t) and u 2 (t), 
since in this case b(t) = e~ 2 4 = W 1 ^ 2 and thus the change of variables in time in (44) 
is trivial, t' — t (as explained in Sec. 2.4). Thus, for the case of the damped harmonic 
oscillator the Generalized Dodonov-Man’ko invariant Xqdm coincide with Xdm- This is 
a special feature of the damped harmonic oscillator, since / 2 ( 0 ) = 0 in this case. 

For an arbitrary wave function ip ho of the harmonic oscillator, the corresponding 
wave function Pck in the Caldirola-Kanai system is given by: 

1 

Vb 


. , . ^ ■ ro b 2 

( pcK{x,t ) = — -j=.e 2h wb 


< Pho 


X X 2 : t -i—'ye^x 2 

t = e 4 e *« 7e x ip HO 
b ' 


2 + 
e 2 x,t 


(50) 


This makes apparent the general strategy: make computations in a (simpler) system 
(e.g. compute p>ho) and map them to the system of interest (get Pck)- 

It should be stressed that in deriving (50) it has not been necessary to solve an 
eigenvalue equation for an invariant in the Generalized Caldirola-Kanai variables, (50) 
is the result of a unitary transformation between two Hilbert spaces. This is a great 
improvement with respect to the Lewis-Riesenfeld or the Dodonov-Man’ko invariant 
method, where first an invariant is found and then its eigenvectors are computed in 
order to provide a basis of the Hilbert space. With the QAEPT it is possible to map 
any harmonic oscillator wave function ipuo into its corresponding function <pcK hi the 
Caldirola-Kanai Hilbert space. In particular number states, coherent states, squeezed 
states, or even density matrices (see [25]). 


3.2. Hermite oscillator 

Similar considerations can be made in more general systems following the same steps. 
Let us consider a LSODE system that has a damping rate linear in time 7 = at, with 
a > 0 : 


x + atx + u> 2 x = 0 . 


(51) 
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This equation is similar to the Hermite differential equation, thus this system is known 
as Hermite oscillator [49]. In analogy with the oscillator with time-dependent mass, the 
Hermite oscillator would have a mass 

m(t) = me“‘^ 2 . (52) 


In order to seek for wavefunctions of the quantum Hermite oscillator we need the 
solutions of the classical equation (51). Those satisfying the canonicity conditions (10) 
are: 

ui(f) = 1 1 F 1 Q + -at 2 / 2 ) , u 2 (t) = iZ\ i; -at 2 /2^) (53) 

and the Wronskian is W(t) = e~ at2 1 2 . Here \F\{a\b\z) denotes the confluent 

hypergeometric function [39]. 

From these solutions the Arnold transformation (7) is given by: 


r = 


m(t) _ 1 lFl (l + fed; -at 2 / 2 ) 

W) = lU (£; 2 ; -oi 2 /2) 

X X 


(54) 


K = 


u 2(t) iFi(g;i;-at 2 /2) ’ 

and with this the QAT (12) is easily derived. 

The solution u 2 {t) has zeros if a < co 2 (underdamping) and has no zeros if 
a > lo 2 (overdamping) or a = u 2 (critical damping). In this last case u 2 (t) = e~^ t2u}2 


and u\{t) = 


Erfi 


Ut \ 


j. Thus, for overdamping and critical damping the 
Arnolnd transformation is defined for all times, whereas for the underdamping case the 
transformation is local in time, mapping a patch of the Hermite oscillator onto a patch 
of a free particle trajectory. Note that in the underdamping case the system performs a 
finite number of oscillations, since the zeros of the confluent hypergeometric functions 
are finite: u 2 (t) is even and has 2f ^~ a ] zeros whereas u\{t) is odd and has 2\ u>2 ~ 2ot ~\ +1 
zeros. 

The Arnold-Ermakov-Pinney transformation (and its corresponding quantum 
version QAEPT) is also easily derived, resulting in: 


t! = 


x = 


W{t) 

w 


dt — — arctan ojqt 
(U 0 


X 


b(t) 


(55) 

(56) 


where cuq is the (arbitrary) frequency of the auxiliary harmonic oscillator, and 


b(t) = u 2 (t)-\Jl + uJqT 2 = \Ju 2 {t) 2 +u^u 1 {t) 2 (57) 

Note that, as before, b{t) never vanishes. 

With this the QAEPT is obtained and the generalized Lewis invariant Zl (i.e. the 
Hamiltonian for the harmonic oscillator mapped into the Hermite oscillator) can be 
computed, recovering the results in [49]. Let us give, instead, the construction of a 
generalized Dodonov-Man’ko invariant as proposed in Sec. 2.3 and 2.4. 
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Choosing h(t) = W(t) = e “* 2//2 and renaming x' 


into 


y + (^ 2 - \ - \^ 2 t 2 )y = o, 


and the wave functions transform as given in (36): 


0 (x, t) = e* at e 


at Z -—i^rrOite 2 X Z 


V(e^ t2 x,t), 


y, the equation (51) transforms 

(58) 

(59) 


where <p(y,t) is a solution of the Schrodinger equation for the time-dependent harmonic 
oscillator with frequency uj 2 — ja 2 t 2 . Applying the Lewis-Riesenfeld invariant 

method to this time-dependent harmonic oscillator we obtain a Lewis invariant that, 
when mapped to the original system provides a generalized Dodonov-Man’ko invariant, 
given by: 

2Adm = 2^P 2 + 2 m ^ uj2 ~ a )X 2 (60) 

Note that there is no term in xp + px ; due to the fact that in this case /(0) = 0. 
This term would appear if we had chosen another initial time i 0 p 0. 

A similar construction can be performed for the Lane-Endem oscillator: 

x-\ - ——x + UnX = 0 , (61) 

1 + vt 0 ’ v ’ 


for which the construction of the Lewis invariant was given in [50]. We shall only provide 
here the expression of the Generalized Dodonov-Man’ko invariant, which turns to be: 


Tgdm 


1 

2m 


P 2 + ^m(oj 2 + 1 yiv)X 2 + 


(i XP + PX 
2 2 


(62) 


4. Conclusion and outlook 


In this paper the relation between the Lewis-Riesenfneld invariant method and the 
Quantum Arnold-Ermakov-Pinney transformation has been established. The former 
aims at finding an invariant for a time-dependent harmonic oscillator, in order to build 
a basis of eigenstates of this invariant satisfying the Schrodinger equation (up to a time- 
dependent phase). The latter is a unitary transformation that relates two Generalized 
Caldirola-Kanai systems, and allows to map states and invariant operators from one 
system to the other. In particular, if one of the system is a harmonic oscillator and the 
other is a time-dependent harmonic oscillator, the Lewis-Riesenfeld invariant method 
is recovered. The time-dependent phase is built into the transformation, in such a way 
that it maps solutions of the Schrodinger equation for one system into solutions for the 
Schrodinger equation of the other directly. 

Any invariant of the harmonic oscillator is mapped, through the QAEPT, into an 
invariant of the time-dependent harmonic oscillator. In particular, the Hamiltonian for 
the harmonic oscillator is mapped to the Lewis invariant, explaining why it has discrete 
eigenvalues. 
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In this paper we have also shown that the QAEPT is global, in contrast to the 
QAT, which is local in time. This explains the robustness and wide applicability of the 
Lewis-Riesenfcld invariant method. 

The main advantage of the QAEPT is that it can also be applied to damped systems 
(or with time-dependent mass) in a unified way. In the case of constant damping 
and frequency we have shown that a more convenient choice is the Dodonov-Man’ko 
invariant, rather than the usual Lewis invariant. The reason is that it shares with the 
Hamiltonian the loss-energy property [18] (which amounts to the only explicit time- 
dependence being through the Wronskian W(t)), and therefore its eigenstates or the 
associated coherent states have nicer properties||. 

We thus propose in the general case to use, as an alternative to the Lewis invariant, 
the generalized version of the Dodonov-Man’ko invariant when dealing with time- 
dependent damped (Generalized Caldirola-Kanai) systems. This could be useful in 
designing short-cuts to adiabaticity processes [16] when damping or time-dependent 
masses are present [35, 34], as well as in constructing coherent and related states. 

Another interesting application of the QAEPT would be in the case of mixed states 
(see [25] where the QAT was discussed in this setting). Let us denote by pno a density 
matrix for the harmonic oscillator, satisfying the quantum Liouville equation 

d -^ = - 1 l [H„o,ho\. (63) 

Then, the density matrix pgck = EpnoE~ l is a proper density matrix (since E is 
unitary Tt(pgck) = Tr (pho) = 1 ) satisfying the quantum Liouville equation for the 
Generalized Caldirola-Kanai oscillator: 

^ = E{Hock,Pgck]. (64) 

All properties of the density matrix puo are transferred to pgck , such as 
characteristic functions, quasi-probability distributions, etc. Also, since Tr (pgck) = 
Tr(p^o), the purity or mixed-state character of p HO is shared by Pgck- fa particular, 
if pho describes a Gaussian state, pgck also corresponds to a Gaussian state. 

However, since the QAEPT does not transform Hamiltonians into each other (only 
Schrodinger and quantum Liouville equations do), care should be taken in the physical 
interpretation of the transformed density matrix. For instance, a thermal equilibrium 
state for the harmonic oscillator is not mapped into a thermal equilibrium state of the 
Generalized Caldirola-Kanai oscillator (in the case when it makes physical sense, for 
instance when the time scale of the time-dependence of the Hamiltonian Hqck is much 
larger than that of relaxation to thermal equilibrium). 

A deeper study of the QAEPT applied to mixed states in order to analyze how 
entanglement is transformed under E (generalized to multipartite systems), and how 
it can be used to describe dissipation and decoherence analyzing the transformation 

|| There is an interesting connection between the Dodonov-Man’ko invariant and the Hamiltonian for 
the damped harmonic oscillator in expanding coordinates H exp , see [47, 48, 51]. 
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properties of master equations like the Lindblad one under E is the subject of a work 
in progress and will be presented elsewhere. 
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